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The use of highly charged ions increases the precision and resolving power, in particular for short- 
lived species produced at on-line radio-isotope beam facilities, achievable with Penning trap mass 
spectrometers. This increase in resolving power provides a new and unique access to resolving low- 
lying long-lived (T1/2 > 50 ms) nuclear isomers. Recently, the 111.19(22) keV (determined from 
7-ray spectroscopy) isomeric state in '^^Rb has been resolved from the ground state, in a charge state 
of g = 8-1- with the TITAN Penning trap at the TRIUMF-ISAC facility. The excitation energy of the 
isomer was measured to be 108.7(6.4) keV above the ground state. The extracted masses for both 
the ground and isomeric states, and their difference, agree with the AME2003 and Nuclear Data 
Sheet values. This proof of principle measurement demonstrates the feasibility of using Penning 
trap mass spectrometers coupled to charge breeders to study nuclear isomers and opens a new route 
for isomer searches. 

PACS numbers: 21.10.Dr, 21.10.-k, 82.80.Qx, 82.80.Rt 



I. INTRODUCTION 

Many advances in our understanding of nuclear struc- 
ture have come from studying nuclei at the limits of 
existence, such as near the particle drip-lines. In very 
neutron- rich nuclei the discovery of nuclear halos 
neutron-skins 0, Q , and the emergence of new sub-shell 
closures and magic numbers Q have provided stringent 
tests for nuclear models. The emergence of new magic 
numbers gives rise to the possibility of low-lying long- 
lived (T'1/2 > 50 ms) isomers, such as spin trap and K 
isomers [5], in neutron-rich nuclei @. Nuclear isomers 
are of interest for a variety of reasons. In astrophysics 
isomers play a role in determining the abundances of the 
elements, since low lying isomers can be excited through 
thermal excitations in hot astrophysical environments 7J , 
e.g. the abundance of ^^Al Q. Furthermore, isomers 
near new magic numbers in neutron-rich nuclei may play 
an important role in determining the r-process path. Iso- 
mers also provide sensitive tests for nuclear models. The 
life-time of a nuclear state depends on the overlap be- 
tween the excited and ground states, a small difference 
in a model wavefunction can lead to widely different pre- 
dictions for nuclear half- lives |9j. Additionally, nuclear 
isomers are intimately tied to nuclear structure through 
nuclear shapes and high spin states Q . 
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The production rates of very neutron-rich nuclei tend 
to be very low, and the possibility of half-lives greater 
than « 50 ms limits the potential to study these iso- 
mers in traditional /? and 7-ray spectroscopy experi- 
ments. Another hindrance for traditional spectroscopy 
experiments is the potential for contaminants or large 
backgrounds ^lOj from molecular beams of lighter nuclei. 
A viable alternative to study the excitation energies of 
these isomers, with lifetimes in the millisecond regime, 
is through mass measurements. Two mass measurement 
techniques that can be used to measure nuclear isomers 
are storage rings and Penning traps [llj. The ability to 
study nuclear isomers has been demonstrated at the GSI 
storage ring, with the discovery of long-lived isomers in 
neutron rich hafnium and tantalum isotopes 12| and in 
the proton rich nucleus ^^^Ce [H, [i3|. Storage rings can 
measure isomer excitation energies as low as 100 keV flU, 
however, the technique is limited to isomers with half- 
lives greater than several seconds due to the length of 
the cooling process. 

Penning traps have been shown to be the most pre- 
cise mass spectrometers for stable 15 1 and unstable 
isotopes. This property along with the ability to per- 
form measurements with low count rates, as low as a 
few ions per hour as demonstrated by SHIPTRAP [IJ, 
allows nuclear structure to be studied in nuclides and iso- 
mers near the particle driplines. To facilitate the study 
of isomers in Penning traps, it is desirable that the time 
of flight resonances, obtained by the time-of-flight ion- 
cyclotron resonance technique |18j , be separated by more 
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than one full-width half-maximum (FWHM) of the reso- 
nance line shape. The FWHM /S.upwHM is proportional 
to the inverse of the radio-frequency excitation time T^f 
as lS.vpwH]\i oc T^p [T^. The mass of the ion is obtained 
from the cyclotron frequency Vc which is extracted from a 
fit of the theoretical line shape. The cyclotron frequency 
of an ion in a homogeneous magnetic field is related to 
its mass as, 



where q and mio„ are the charge and mass of the ion 
of interest and B is the magnetic field strength in the 
trap. Thus, the resolving power TZ required to separate 
an isomer from the ground state by one FHWM is. 

It — — oc — (2) 

ArUa ZTTnia 

where rUa is the atomic mass of the ground state of the 
nuclide of interest and Ama is the difference in atomic 
masses of the isomer and ground state. In what follows 
resolved will be taken to mean separated by one FHWM 
of the line shape. Since the charge state and the excita- 
tion time in the Penning trap both enter into Eq. ([2]) in 
the numerator a decrease in the excitation time can be 
compensated by the appropriate increase in the charge 
state. As a rule of thumb the excitation time is chosen 
according to the relation T^p < 3 ■Ti/2- 

The ability to resolve isomers in Penning traps was first 
demonstrated at ISOLTRAP, where the isomeric states 
of *^'^*Rb were seen [l9|. If the hyperfine structure of a 
nucleus is known, then through state selective laser ion- 
ization it is possible to create isomerically pure beams 
from which the ground and isomer masses can be deter- 
mined with a minimum of contamination (20l . [2l| . If the 
hyperfine structure of the ground and isomeric states lie 
close to each other then it is possible to create pure beams 
by the use of an appropriate ion/Penning trap cleaning 
technique [13, H^l ■ The shortest of these cleaning tech- 
niques still requires upwards of one-hundred milliseconds 
and cannot be used for nuclei and isomers with half-lives 
on the order of tens of milliseconds. Penning traps can 
also be used as a discovery machine with the first ob- 
servation of new nuclides [2^ or isomeric states (13], or 
to determine isomer excitation levels seen in decay spec- 
troscopy im,!!!,!!!]. 

The increase in resolving power from charge breeding 
enables the resolution of low lying nuclear isomers in Pen- 
ning traps. The dots in Fig. [T] show the normalized rela- 
tive difference in mass Am/m? for all known isomers with 
half-lives greater than 1 /is and excitation energies less 
than 9 MeV (868 in total) plotted against the half-life of 
the isomer. Superimposed as lines are calculations of the 
available mass normalized resolving power for a given q 
and charge breeding time Tft^eed assuming an excitation 
time of 3 • — Tbreed- The hatched area highlights the 
gain in resolving power, as these isomers could not be re- 
solved as singly charged ions due to their short life-times. 




FIG. 1. (color online) The normalized relative mass difference 
between an isomer and the ground state plotted against T-1/2 
of the isomer. Superimposed is the available mass normalized 
resolving power, Eq.J^ for different 5's plotted against Trf- 
Data are from Ref. [27l |. The vertical line at 3 ms represents 
the shortest half-life that may be measured at TITAN. See 
text for details. 

Isomers with half-lives much shorter than the breeding 
time required to reach a given charge state cannot be 
resolved. Any isomer to the right of the lines can, in 
principle, be resolved provided sufficient yield and a rea- 
sonable isomer to ground state production ratio. Ideally, 
in the absence of an isomer purification mechanism, both 
states would be delivered in equal amounts. A ratio con- 
siderably different than this would cause the resonance of 
the lesser produced species to be lost in the background 
produced by the primary species. In Fig. [1] isomers to 
the left of a line cannot be resolved with the chosen q. 
To demonstrate the improvement in the resolving power 
gained by using highly- charged ions to study nuclear iso- 
mers a mass measurement of ''^Rb was performed in a 
charge state of g = 8-1- and the reference ion was *'^Rb^+. 
For the first time the ground and isomeric states of ^^Rb 
were resolved such that the minima were separated by 
more than one full-width half-maximum. 



II. THE TITAN FACILITY AND 
EXPERIMENTAL SET-UP 

Triumf 's Ion Trap for Atomic and Nuclear physics (TI- 
TAN) [iiliil is located in the IS AC haU of the TRIUMF 
laboratory. The TRIUMF-ISAC facility [30,] produced 
the ^^Rb beam by impinging a 500 MeV proton beam 
with a current of up to 98 /iA beam on a Nb production 
target. The produced radionuclides thermally diffuse out 
of the target where they are ionized and accelerated to an 
energy of 20 keV. The beam is cleaned of contaminants 
using a 120° dipole magnet with a resolving power of 
m/5m -3000 and delivered to TITAN. TITAN currently 
consists of three ion traps: a radio-frequency quadrupole 
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beam to next experiment 



off-line ion source 



FIG. 2. (color online) Illustration of the TITAN experimental 
set-up. For precision mass measurements singly charged ions 
delivered from ISAC or the offline source can (a) be injected 
directly into the MPET or (b) first sent to the EBIT for charge 
breeding and then sent to the MPET. 



(RFQ) cooler and buncher [3i|, an electron beam ion 
trap (EBIT) 13211 and a high precision measurement trap 
(MPET) QlHil. A schematic of TITAN is shown in Fig. 

El 

The RFQ is used to stop, thermalize, and bunch the 
beam delivered from ISAC in a helium buffer gas. This 
cooling decreases the energy spread of the beam which 
allows for a more efficient injection into either the MPET 
or EBIT. The beam is extracted from the RFQ with an 
energy of 2 keV. The cooling is of chief concern for the 
MPET since a large energy spread adversely affects the 
precision obtained in a mass measurement. 

The MPET accepts beam either directly from the RFQ 
as singly charged ions (Fig. [2l path (a)) or from the 
EBIT as highly charged ions (Fig. [21 path (b)), and 
performs high precision mass measurements by utilizing 
the timc-of-flight ion cyclotron resonance technique (see 
Ref. [l8|). A magnetic field of 3.7 T in the MPET is 
used to radially confine the ions while electrostatic fields 
are used to trap the ions axially. The general procedure 
for mass measurements with TITAN is outlined in [33j . 
The TITAN system is a well established measurement 
facility for singly charged ions, in particular for very- 
short lived, iieutron-rich light isotopes such as ^'®He [3^, 
Hi, "Li [33, and "'i^Bg g Recently, the first 

mass measurements of short-lived, highly charged ions 
were carried out at TITAN 

The TITAN EBIT charge breeds the ions of interest 
stored initially as singly charged ions in the trap through 
successive impact ionization of the electron beam. The 
charge breeding process increases the ions' energy spread, 
which is estimated to be in the range of 10 — lOO-g eV 41|. 
In order to cool the charge bred beam a cooler Penning 
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FIG. 3. (color online) Time-of-flight spectra of charge bred 
ions from the EBIT detected on a micro-channel plate before 
the MPET. The red (black) curve shows the case of injection 
(no injection) of *^Rb into the EBIT. 



trap is currently being built and will use either electron or 
proton cooling to decrease the phase space of the highly 
charged ion bunch [42, 4^] . A full description of the TI- 
TAN EBIT can be found in Ref. [13, Hi. The ions are 
radially bound by a 3 T magnetic field and the space 
charge from the electron beam, and confined axially by 
a 100 V deep well. The injected ions and atoms of resid- 
ual gas in the trapping region are bred to higher charge 
states and are extracted from the trap with an energy 
of E]^^^ « 1.9 • q keV. The EBIT is typically opened for 
a short time, between a few hundred nanoseconds to a 
few microseconds, allowing for the selection of ions of 
different m/q' s by time-of-flight gating using a Bradbury 
Nielson gate [4^. For this experiment the ^*Rb beam 
was charge bred for 23 ms with an electron beam energy, 
relative to the trapped ions, of 2.5 keV and a current of 
10 mA, yielding a similar charge state distribution pre- 
sented in Fig. [3] for ^'^Rb. An m/q ratio of ~ 9.5, which 
corresponds to ionic states of ^^Rb^"*" and ^^Rb®"*", was 
chosen due the cleanliness (lack of charge bred residual 
gas) of the TOF spectrum, hence, a clear separation of 
species could be achieved. 



III. DATA ANALYSIS 

The quantity of interest in a Penning trap mass mea- 
surement is the cyclotron frequency i^c of the trapped 
ions, and is related to the mass as shown in Eq. ([T]). This 
is extracted from a fit of the theoretical line shape [3| 
to the time-of-flight resonance spectra. To eliminate ma- 
jor systematic shifts, such as those due to drifts of the 
magnetic field, the ratio R of the cyclotron frequency 
of the ion of interest to that of a well known reference 
ion is taken in the combination such that the reference 
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FIG. 4. Illustration of the correlation introduced between ad- 
jacent frequency ratio measurements due to shared references. 
The filled circles (tj's) represent reference measurements of 
Vc,ref and the open circles (Ti's) show the interpolation of 
Vc,ref to the center time of a measurement of an ion of inter- 
est. From the figure it is clear that Ti is correlated with Ti+i, 
Ti and Ti+i with Ti+2, and Ti+2 with Ti+s (details follow in 
the text). 



frequency is in the numerator, 



'-^c.ref 



(3) 



Taking the ratio in this way greatly simplifies the calcula- 
tion of the covariance relations presented in Sec. (IIII Al) . 
After averaging a number of frequency ratio measure- 
ments Eq. (ID) can be used to determine the mass of the 
ion of interest relative to the mass of the reference ion, 



rUa = R{ma,ref 

Qref 



(4) 

where R is the weighted average of all frequency ratios 
and Be and B^^ref are the electron binding energies of 
the ion of interest and of the reference ion. The electron 
binding energies for Rb^+^^+ are w 500 eV and 650 eV 
p6j . The population of long-lived ionic metastable states 
that add significant amounts of energy, or equivalently, 
mass, to the system can be neglected, as the binding 
energies are relatively small. 



for isomers or nuchdes that lie very close in mass and are 
trapped quasi-simultaneously, as Coulomb interactions 
between the two stored species in the trap lead to system- 
atic shifts in the measured frequencies [l^, 113] . Therefore 
a linear interpolation between two frequency measure- 
ments, one taken immediately before and one immedi- 
ately after the ion of interest must be made. As a time 
saving measure during an experiment two measurements 
of the ion of interest will often share a reference mea- 
surement introducing correlations between the frequency 
ratios. With the use of highly charged ions and the high 
level of precision that can be reached it is important to in- 
clude these correlations when determining the final aver- 
aged frequency ratio. The relative statistical uncertainty 
of the cyclotron frequency in a measurement is related to 

the resolving power as Sma/ma oc TZ~^ • ^/N [i^ where 
N is the number of detected ions. Here we present two 
cases shown in Fig. |4l First, the most likely case where 
two measurements share a reference measurement, and 
second, the case where several measurements occur be- 
tween two reference measurements. In practice the sec- 
ond case does not occur since these data are generally 
summed. In the second case, the analysis with the time 
correlations and with the summed data will yield nearly 
the same result because the summed data implicitly in- 
cludes time correlations between the frequency measure- 
ments. For the first case, the covariance relation between 
frequency ratios is. 



covar(i?.j,i?.j+i) 



1 



T- -t- 



i+2 



(5) 

where the i and i + 1 refer to the ith and «th+l mea- 
surements of the ion of interest and the j's refer to the 
reference measurements. For the second case, the covari- 
ance between frequency ratios that share both references 
is. 



covar(i?j,i?j+i) 



1 



J+1 



T- 



1 - tj 



T- 



i+1 



(6) 



Figure|4]illustrates the relationship between the variables 
given in the above equations. In both cases the covari- 
ance is proportional to the variance of the reference mea- 
surements. It is desirable to measure the reference ion 
much more precisely than the ion of interest to reduce 
correlation effects, however, a trade-off must be made in 
order to maximize the statistics collected, and hence, the 
precision, for the ion of interest. 



B. Systematic errors/Uncertainties 



A. Time correlations between data sets 

It is generally not practical to measure the reference 
frequency at the same time as the ion of interest, except 



Several systematics must be taken into account. Sys- 
tematics relating to misalignment between magnetic and 
trap axes, electric field miscompensation, relativistic ef- 
fects, etc., are minimized by choosing a reference ion 
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FIG. 5. (color online) Time of flight spectrum of ions ex- 
tracted from the MPET. See text for details. 



which is close in m/q to the ion of interest as these ef- 
fects scale with the difference in the charge to mass ratio 
A{m/q) [ll]. To determine any potential shifts due to 
different m/q effects between the ion of interest and the 
reference ion, a series of mass measurements on 85j^|^io,8+ 
and ^^Rb^"*" using ^^Rb®"*" as the reference were com- 
pleted. The extracted masses all agree within 1 cr of the 
literature value. Although no shifts were observed to be 
conservative we take, as an upper limit on any systematic 
effects, a systematic uncertainty of 42 parts-per-billion 
(ppb) in the frequency ratio. 

A second systematic effect stems from the ambiguity 
in selecting the upper and lower time cuts on the time of 
flight spectrum. The ambiguity arises due to charge ex- 
change processes in the trap. If an ion undergoes charge 
exchange with residual gas in the vacuum, these ions will 
manifest themselves as a long tail in the time of flight 
spectrum. In the present analysis the lower and upper 
levels were set at 12 and 40 /is, respectively. The lower 
level was set to 12 /is in order to maximize the num- 
ber of on resonance ions while minimizing background 
counts from the nearby peak resulting from charge 
exchange in the trapping region. Figure [5] shows a typical 
time of flight spectrum for ^^Rb^"*" which was trapped for 
197 ms. The dashed-blue lines show the lower and upper 
time cuts while the solid-red lines show when, on aver- 
age, ^*Rb ions with different charge states would arrive. 
In order to determine the systematic effect R was de- 
termined for upper level time cuts of 30, 35, 40, 45 and 
55 /zs for both the ground and isomeric states. If the av- 
erage frequency ratio determined at 40 /Lts for either case 
was an extremum the systematic effect was assigned to 
be the full range of the extracted ^'s, otherwise half of 
the range was assigned. 

The last systematic effect results from ions with dif- 
ferent m/g's which are in the trap at the same time. 
For ions that are separated by more than the full-width 
half-maximum of the resonance curve the measured cy- 



clotron frequency of both ions decreases from the nomi- 
nal value [l9[. In order to eliminate this potential shift 
a count class analysis [13] was performed. For this the 
frequencies are determined as a function of the num- 
ber of detected ions, hence the number of ions stored 
in the trap (count class). The cyclotron frequencies of 
the ground and isomeric states were extracted for differ- 
ent count classes and the frequencies were extrapolated 
to the detector efficiency of 0.6 ± 0.2%. When one ion 
was detected after extraction on average about 1.7 ions 
were actually in the trap, thus, an extrapolation past 
unity was required. The count class analysis eliminated 
any correlations between the extracted fit parameters, so 
their correlations were not included in the analysis. 

C. Extracting isomeric excitation energies 

The excitation energy of the isomer can be extracted 
by using Eq. and, using the isomer as the reference, 
the difference in mass is 

rUaAso - ma,gnd = (1 " R){ma,iso ~ qrUe + Be). (7) 

Since the ground state and isomer represent a mass dou- 
blet, many of the systematic effects presented in the pre- 
vious section cancel. However, to be conservative we in- 
clude the maximum 42 ppb systematic shift presented in 
the previous section. 

If other reference measurements were performed, it is 
also possible to extract the mass excess of each state as 
well as the energy difference of the isomeric state through 
a joint fit of all the frequency ratio pairs of the ground 
and isomeric states, making use of the full covariance 
matrix between all the frequency ratios using Eqs. [5] 
and [6l This process provides a better error estimate on 
the difference than achieved using simple error propaga- 
tion methods. The weighted average of each frequency 
ratio, to be used in Eq. IH can be extracted from the 
same formalism and will yield a higher uncertainty than 
a weighted average assuming independent measurements 
since the correlations are positive. The technique used is 
fully described in Ref. 

IV. RESULTS 

The frequency ratios for ^®Rb®+ (ground 
state) and 7*^"Rb8+ (isomeiB) (T1/2 
17.66(3) min, 5.74(3) min [5T|) measured relative 
to *^Rb^+ are presented in Table [D The error budget, 
showing the contributions from all systematics discussed 



We use the notation introduced in the NUBASE |5C| evaluation 
where isomeric states are labeled by m, n, p, q, etc. in order of 
increasing excitation energy. 
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TABLE I. Frequency ratios and mass excesses of "8. '»"]^i3»+ measured relative to Rb with the total uncertainties. The 
error budget is summarized in Table Ull 



Isotope 




ME (keV) 


MEameoz (keV) 


ME-AIEameos (keV) 


''«Rb«+ 


1.032475265(99) 


-66933.2(7.0) 


-66936.228(7.452) 


3.0 




1.032476806(60) 


-66824.8(4.2) 


-66825.038(7.455) 


-0.2 
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FIG. 6. (color online) Resonance of Rb for quadrupole 
excitation times of (a) 97 ms and (b) 197 ms. The isomer 
resonance is on the left and the ground state resonance is on 
the right. The solid line is a fit of the theoretical line shape 
[l^ to the data. 



TABLE II. Error budget of the frequency ratios in ppb. 
"Ground" and "Isomer" indicate the measurements were done 
relative to *^Rb^^, while "Doublet" indicates that the isomer 
was used as the reference ion. 



Uncertainty (ppb) 


Ground 


Isomer 


Doublet 


Statistical 


66.8 


29.5 


66.7 


Time Correlations 


24.1 


24.5 


0.0 


TOF Spectrum Cut 


55.0 


20.0 


38.0 


Trap Systematics 


42.0 


42.0 


42.0 


Total 


99.0 


60.0 


87.0 



TABLE III. Extracted excitation energy of the isomeric state. 
The error in parentheses is the statistical and time correlation 
error, while the error in braces is the total error. The time 
correlations in this measurement are of the second type (Eq. 

Method Excitation Energy (keV) 

Time Correlation 108.4(4.8){7.6} 
No Time Correlations 108.5(5.1) {7.8} 
78n-p^^8+ ^ reference 108.7(4.8){6.4} 



previously, is presented in Table |TT1 Typical resonances 
for 97 ms and 197 ms excitation times are shown in 
Fig. ini Subtracting the literature excitation energy 
of the isomer from the measured mass excess and 
taking the weighted average of both measurements (at 
Trf = 97, 197 ms), it was found that the mass excess 
for ^^Rb is -66935.3(3.6) keV. This is in excellent 
agreement with the AME03 value [55|. 

Extracting the mass of the ground state was com- 
plicated by two factors. First, the ratio of iso- 
meric to ground state nuclei delivered from ISAC was 
Nisomer/Ng,s. ~ 2. The difference in the number of ions 
weakens the time-of-flight resonance of the ground state 
leading to a larger uncertainty in the fitted frequency. 
Second, while the resonances are fully resolved in Fig. 
[6l[b), the ground state sits very near the first sideband 
of the isomer resonance leading to the ambiguity of what 
the fitting routine is extracting: The ground state cen- 
ter frequency or the position of the first sideband of the 
isomer? In our case the sideband of the isomer is half 
of the depth of the ground state resonance causing the 
isomer sideband to have an affect on the fitted ground 
state frequency. 

For the TITAN set-up full resolution of the isomer and 
ground states is achieved with a quadrupole excitation 
time of 197 ms shown in Fig. [S] (b). The excitation 
energy of the isomer when using ^^"Rb^+ as the refer- 
ence is 108.7(6.4) keV. When a joint fit is performed 
on the data, taking into account the correlations be- 
tween measurements, the energy of the isomer relative 
to the ground state is 108.4(7.6) keV. When time corre- 
lations are neglected, the excitation energy of the isomer 
is 108.7(7.8) keV. The results are summarized in Table 
mil Both are in agreement with the more precise value 
of 111.19(22) keV obtained from 7-ray spectroscopy [5l|. 
It is unsurprising that the excitation energy derived from 
using ^^"Rb*"*" as the reference is more precise since fewer 
systematic effects enter into the calculation. It is interest- 
ing to note that the uncertainty on the energy difference. 
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before systematic effects are included, for both the case 
with time correlations and with ^^"Rb*+ as the reference 
are identical. 

V. CONCLUSION 

The ground and isomeric state in ""^Rb have been fully 
resolved in TITAN's precision mass measurement Pen- 
ning trap through the use of highly charged ions. This is 
the first time a full separation of the resonance shapes, for 
this isotope and isomer, could be achieved in a Penning 
trap. The mass excess values extracted for both states 
agree with the AME2003 values, as does the measured 
excitation energy of the isomer. The uncertainty on the 
weighted ground state mass excess has been reduced by 
a factor of 2 compared to AME03, resulting in a value 



of —66935.3(3.6) keV. This is a proof-of-principle that 
highly charged ions are a powerful tool for increasing the 
precision and resolution of Penning trap mass spectrom- 
eters for on-line spectroscopy and searches for isomeric 
states. 
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